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ABSTRACT 

> 

CNj The trace anomaly is relevant to the renormalisation of the unpolarised parton distributions. The 
coefficient which describes the local 7* (glue) interaction differs at next-to-leading order between the 
renormalised operator product expansion and parton model descriptions of deep inelastic scattering - 
independent of the factorisation scheme. 
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There are two approaches to the theory of deep inelastic scattering (DIS). These are the operator 
product expansion (OPE) and the QCD improved parton model (PM) [1]. In both approaches we write 
the transverse and longitudinal structure functions as the sum over the convolution of soft quark and 
'gluonic parton distributions with hard coefficients, viz. 

F[T,L](x,Q'^) = + ® C'[T,Li(^'«*'Q^) + xg (S)C^j,^^{x,as,Q'^) (1) 

The hard coefficients C^, and C^, describe the interaction of the hard photon (carrying transverse 
or longitudinal polarisation) with the quark and gluon partons. 

In the parton model we calculate the deep inelastic cross section from a quark or gluon target 
to some order in perturbation theory and then use the factorisation theorem to absorb the infra-red 
coUinear singularities into the parton distributions — or, equivalently, the hadron wavefunction. The 
remaining infra-red safe part of the phase space is identified with the hard coefficients. This factorisation 
is universal to all inclusive, high-energy hadronic interactions. 

Having taken care of the collinear singularities we also need to renormalise the parton distributions. 
Here we use the operator product expansion, which is applied to the total deep inelastic cross section. 
The operator product expansion defines the target wavefunction in terms of the matrix elements of a 
tower of renormalised, local operators [1]. These operator matrix elements need to be renormalised 
consistent with the symmetry constraints of anomaly theory [2,3] — in this case, the trace anomaly [4]. 
This renormalisation has to be carried out in addition to the renormalisation of the propagators and 
vertices in QCD. In this paper I explain how the anomaly induces a difference at next-to-leading order 
between the quark distributions and gluonic coefficients which are obtained via the operator product 
expansion and the (factorisation only) parton model. This difference is independent of the factorisation 



itself in the gluonic coefficient. Finally, these results are compared with the axial anomaly in polarised 
deep inelastic scattering. 

The parton distributions in equ.(l) arc defined with respect to the operator product expansion so 
that their even moments project out the target matrix elements of the renormalised, spin-even, local 
operators 



= < p\ 

k 



R 



\P >c (2) 



and 



ip+y^+' f dx x'^+'gix, Q') =< p\ [lVG+«(0)(zD+)2"G\(0) 
Jo L 



\P >c (3) 



Here G^, is the gluon field tensor, = + igA^ is the gauge covariant derivative in QCD and is 
the proton momentum vector. The second (n = 0) moment of the quark (gluon) distribution projects 
out the target matrix element of the quark (gluon) part of the energy- momentum tensor O^j,^ in QCD. 
The operators in equs.(2,3) should be quoted with respect to some renormalisation prescription it! — 
here we have written them at the subtraction point /x^ = Q^- 

The problem of anomalies is related to the symmetry of the Dirac sea at infinite momentum. When 
we quantise the theory about some background mean field (the non-perturbative vacuum) then we 
construct the Fock expansion by keeping quark and gluon quanta with momentum less than some ultra- 
violet cut-off A. High momentum (short distance) fiuctuations greater than this cut-off are absorbed 
into the theory as part of the definition of the mean field. The anomaly has an analogy in classical 
physics where we consider a thermodynamic system in equilibrium with a heat bath [5] . There can be a 
nett fiow of energy between the system and the heat bath at equilibrium so that the hare energy of the 
system is not conserved. However, we can write down a free energy, which is conserved. We may think 
of the perturbative Fock space as the system and the mean field as the heat bath. One finds that there 
is a fiux of energy between the mean field and the Fock space but the free energy remains constant [5] . 
The physical energy momentum tensor corresponds to the free energy. It is conserved but picks up 
a finite trace due to the anomaly [6]: 

= (1 + 7m) E mm + ^G,^G>^r (4) 

Here ruq is the running quark mass, 7^ is the mass anomalous dimension and I3{as) is the QCD beta 
function. We can regard the Fock space as a closed system (independent of the mean field) only when we 

consider anomaly free quantities such as the vector current, which measures the valence quark number. 
The anomalous gluonic term in equ.(4) is particularly important in hadronic physics because of the 
contribution it makes to the mass of the proton [7,8]. The matrix element of O^i, in the proton is 

<p\0,j,v\p>^'^P^iPv (5) 

The trace of this equation is 

<p\e^;,\p>=2ml (6) 

If we were to lose the anomaly, then we would break Poincare invariance at the point of regularisation 
and our parton distributions would correspond to a proton with the wrong mass. 

To see how the anomaly contributes to the unpolarised structure function we apply the operator 
product expansion analysis to unpolarised photon gluon fusion at 0{as) — Fig. 1. Here = {qq; Ot, Qs) 
and Pfj^ = {po',0t,P3) are the momenta of the hard photon and the soft gluon respectively; k is the 
momentum of the struck quark or anti-quark. We consider the transverse cross section where the 
photon and gluon both carry transverse polarisation and follow the discussion of polarised deep inelastic 
scattering given in the papers of ref.[9]. A finite quark mass m and gluon virtuality — is used to 



The cross section for this photon gluon fusion process is: 
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where VF^ is the centre of mass energy squared. There are three distinct contributions to the total phase 
space. The first term in equ.(7) receives contributions only from ~ Q^, where the photon makes a 
local interaction with the target gluon. The second term receives contributions from the full range of 
k^. It is sensitive to the mass singularity at zero k'^ and gives the In dependence of the total cross 
section after we integrate over the transverse momentum. The third term is infra-red dominated and 
receives contributions only from k^ ~ m^, —p^. 

At 0{as) the cross section for photon gluon fusion is the sum of the gluonic coefficient C^(x, q;^) 
and the quark distribution of the gluon {q + qY^\x,Q'^). In the parton model we set the coefficients 
equal to the "hard" part of the phase space, where "hard" is defined to include the whole phase space 
which is not subtracted into the target wavcfunction via the factorisation theorem. The local 7* (glue) 
interaction is clearly "hard" in this picture and is included into the parton model as a part of the gluonic 
coefficient C^{x,as). The simplest way to isolate the mass singularity involves choosing a cut-off on 
the transverse momentum of the quarks so that we identify the gluonic coefficient with the phase space 
corresponding to quark jets with transverse momentum k^ > fJ-'jact^ where ^j^ct some finite scale 
(Q^ ^ l-tjact l^}act ^ —p^jfin?). With this cut-off, the hard parton model coefficient is: 
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The remaining soft part of the phase space is identified with the quark distribution of the gluon x{q + 
9)pm(^5 A*/act)' which includes all quarks and anti-quarks with transverse momentum less than /J^'jact- 
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We now compare this parton model distribution with the renormalised quark distribution, which we 
calculate via the operator product expansion. 

Working in the light cone gauge (^4+ = 0), the quark distribution is the renormalised version of the 
graph in Fig. 2: 
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The k^ and k_ integrals are carried out using the two delta function constraints to yield an ultra-violet 
divergent integral in the Euclidean A;|^, which needs to be renormalised. We can choose to regularise this 
ultra-violet divergence via dimensional regularisation or the introduction of a Pauli-Villars term - both 
of which respect the anomaly. If we use dimensional regularisation, then we continue the physics in the 
transverse dimensions from D = 2 to D = (2 — 2e) dimensions; the ultra-violet divergence appears as a 
term proportional to -, which is subtracted out as the infinite renormalisation constant of the parton 
distribution. The sum over the transverse polarisations of the gluon is 



Equ.(lO) evaluates as 
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in D = (2 — 2e) dimensions. If no ultraviolet cutoff is imposed, then the integral 



J k^ + m^ — p^x{l - X) 

develops an ultra-violet pole -, which comes from k"^ much greater than the parton model cut-off 
/Jjacf This pole cancels with the e term in the numerator of equ.(12) to yield a finite term which is 
induced by the regularisation: the (trace) anomaly in the energy-momentum tensor. This anomalous 
term is independent of how we next make the minimal subtraction to obtain the renormalised parton 
distribution. 

In the MS scheme the renormalised quark distribution is 

(q + q)'i\x,fi') = (q + q)%(x,fi') + C'^(x,a,) (14) 
where is the anomaly at 0{as) 

(X 
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TT 

and \j? is the rcnormalisation scale. Here we have sTibtracted out the divergent term ^(^x^ — 2x -\- 
~ 1e + ln47r) into the rcnormalisation constant. In the general minimal subtraction scheme we 
subtract out a divergent term proportional to (^-|-a), where a is some constant; the renormalised quark 
distribution is equal to the parton model distribution in equ.(9) with ii^fact = ^2gin47r-7B-a p^^g -^.j^g 
subtraction independent anomaly, which is left behind by the ultra-violet regularisation. Pauli-Villars 
regularisation yields the same renormalised quark distribution as equs.(14,15). In this case, the anomaly 
comes from k'^ of order the mass squared of the Pauli-Villars fermion in the second term of equ.(12). 

The renormalised quark distribution corresponds to a gluonic coefficient which is different to the 
parton model by the anomalous term in equ.(15); viz. 

g2 Q2 
Cope{^j ~~2 — ' '-^«) ~ ^pm{^j ~~2 — j^s) — C^{x, as) (16) 

t^fact f^fact 

Since the anomaly term C^(x, ag) is a polynomial in x at 0{as), it describes a local 7*(glue) interaction. 
Some of this local interaction is included as a contact interaction with the target wavefunction within 
the OPE approach, which offers the most rigorous treatment of the total deep inelastic cross-section. 
The parton model treatment of the local 7* (glue) interaction is complete if one can treat the proton 
target as an ensemble of perturbative quark and gluon degrees of freedom. However, it is not clear ab 
initio that the parton model "scheme" is sufficiently general that it can describe all local interactions 
between the hard photon and the background colour mean field in addition to the perturbative glue. 
In particular, it is not clear that the former interaction has the same two-quark-jet signature that 
we expect from photon gluon fusion in perturbation theory. The background field via the anomaly is 
important in determining the mass of the proton, equ.(6). We need to understand how it matches onto 
the perturbative glue. The factorisation theorem tells us how to treat the infra-red coUinear singularities 
but does not tell us how to treat the local 7* (glue) interaction, which is strictly non-collinear. 

The parton model gluonic coefficient is recovered from the operator product expansion analysis if 
we define our parton model quark distribution to be 

xiq + q)pMix, Q') = x{q + q)R{x, Q^) - xqr (g) C\{x, Q^) (17) 

fn -1-7?') ^.J^ n'^\ 



is the difference of two renormalised quark distributions — each of which includes the effect of the 
anomaly at infinite momentum. 

In the parton literature it is commonly assumed that the parton model distributions can be written 
as light-cone correlation functions [10,11] 

q(x) = ^ [ dz.e-''''-P+ <p\q{z.)^+Ve''fo~'^y-^+q{0)\p>^ (18) 
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(The moments of these correlation functions project out the target matrix elements of the operators 
in the OPE.) It is clear from our discussion that this identification needs to be modified according to 
equ.(17) because of the anomaly. 

The correlation function offers a simple, physical picture of the anomaly in deep inelastic scattering. 
In a quark model (which has an inherent ultra-violet cut-off) the non-local matrix element in equ.(18) 
has the interpretation that we take out a quark (insert an anti-quark) in the nucleon and re-insert it 
(take it out) at position along the light-cone. The z^ ^ limit of this matrix element is the same 
object that appears in Schwinger's derivation of the anomaly, which uses point splitting [12] (along the 
light-cone). Intuitively, one can think of the anomaly in deep inelastic scattering as a zero correlation 
length effect which is missing in the parton model. 

We now compare our results with the anomaly in polarised deep inelastic scattering. The anomaly 
contribution to the renormalised, unpolarised quark distribution (equ.(14)) arises in the same way as the 
axial anomaly contribution to the spin dependent quark distribution [9], which is measured in polarised 
deep inelastic scattering. The anomaly comes from much greater than the parton model cut-off. 
Clearly, the theory of polarised and unpolarised deep inelastic scattering should be formulated in a 
consistent way. If we imposed some ultra-violet cut-off on the then we would lose the trace anomaly 
in the unpolarised structure function and also the axial anomaly in the polarised structure function. 
This means that the EMC spin effect [13] and the physics of the unpolarised structure function have to 
be discussed together as part of the same problem. 

The renormalised spin dependent quark and gluon distributions which are measured in polarised 
deep inelastic scattering are 

2Ms+(p+)2" f^dx x''''AqJx,Q^) =< - ''^ ^^n^'' 
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where G^j, = \^ij.vapG°'^ is the dual tensor and denotes the proton spin vector (M is the proton 
mass). The spin dependent distributions have to be renormalised consistent the axial anomaly, which 
says that any ultra-violet regularisation that respects gauge invariance does not respect chirality as a 
good symmetry. The renormalisation of each of the C-even axial tensors in equ.(19) involves a gauge 
dependent counterterm A;^^^ ^2 • These fc^^^ ^2 define a gauge dependent gluonic distribution n{x, Q^) 
[14]. 

There is an important physics difference between the anomalies in unpolarised and polarised deep 
inelastic scattering. In the unpolarised case the physical energy-momentum tensor is conserved: the 
canonical identification oi x{q + q)ji with a quark momentum distribution holds true in the renormalised 
theory. The physical axial vector current is not conserved [15] and has a two-loop anomalous dimension, 
which is induced by the anomaly [16] — it does not measure spin [14,17]. This means that AgR(x, Q"^) 
does not have a canonical interpretation. The colour mean field in the proton should be manifest 
explicitly in Aqji{x, Q"^). This effect may be seen as an OZI violation in the spin dependent structure 
function gi at large x [14]. 

Just as in the unpolarised case (equ.l6)), the renormalised quark distribution in equ.(19) yields a 
spin dependent gluonic coefficient which differs from the parton model. The parton model coefficient 



where is the anomalous coefficient at 0{as). It is important to distinguish the Agji{x, Q^) in equ.(21) 
from the gauge dependent k,r{x, Q"^). In the light-cone gauge = (where the parton model is usually 
formulated) one finds kr^x^Q"^) = Ag^i ® C^i^x^Q"^). However, in a covariant gauge kji{x,Q'^) is not 
invariant under gauge transformations — in contrast to the physical Aq^^XjQ'^) and AgR{x,Q'^) [14]. 
In polarised deep inelastic scattering the axial anomaly introduces a second set of (gauge dependent) 
gluonic operators. In unpolarised deep inelastic scattering the trace anomaly involves the same leading 
twist operators that define the regular gluon distribution in equ.(3) in order to yield the trace in equ.(4) 
— see also [5] . 

In summary, the anomaly means that the local 7* (glue) interaction is treated differently in the op- 
erator product expansion and parton model treatments of both unpolarised and polarised deep inelastic 
scattering. Since an operator analysis of anything other than the total deep inelastic cross-section is 
intractable, there is no clear way of renormalising the initial and final state wavefunction in other high- 
energy collisions. It is interesting to consider the possibility that the anomaly might induce a process 
dependence in the parton distributions at greater than or equal to next-to-leading order in a^. 
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